Viscous spreading of an inertial wave beam in a rotating fluid 
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We report experimental measurements of inertial waves generated by an oscillating cylinder in 
a rotating fluid. The two-dimensional wave takes place in a stationary cross-shaped wavepacket. 
Velocity and vorticity fields in a vertical plane normal to the wavemaker are measured by a corotating 
Particule Image Velocimetry system. The viscous spreading of the wave beam and the associated 
decay of the velocity and vorticity envelopes are characterized. They are found in good agreement 
with the similarity solution of a linear viscous theory, derived under a quasi-parallel assumption 
similar to the classical analysis of Thomas and Stevenson [J. Fluid Mech. 54 (3), 495-506 (1972)] 
for internal waves. 



I. INTRODUCTION 

Rotating and stratified fluids both support the prop- 
agation of waves, referred to as inertial and inter- 
nal waves respectively, which share numbers of similar 
properties.!^ These waves are of first importance in the 
dynamics of the ocean and the atmosphere,^ and play a 
key role in the anisotropic energy transfers and in the 
resulting quasi-2D nature of turbulence under strong ro- 
tation and/or stratification.^ 

More specifically, rotation and stratification (and the 
combination of the two) lead to an anisotropic dispersion 
relation in the form a = /(/c^/lkl), where a is the pulsa- 
tion, k is the wave vector, and the z axis is defined either 
by the rotation axis or the gravity.^ This particular form 
implies that a given excitation frequency <j selects a sin- 
gle direction of propagation, whereas the range of excited 
wavelengths is set by boundary conditions or viscous ef- 
fects. A number of well-known properties follow from 
such dispersion relation, such as perpendicular phase ve- 
locity and group velocity, and anomalous reflection on 
solid boundaries 

Most of the laboratory experiments on internal waves 
in stratified fluids have focused on the properties of lo- 
calized wave beams, of characteristic thickness and wave- 
length much smaller than the size of the container, ex- 
cited either from local^ ^^or extended^ ^ sources. On the 
other hand, most of the experiments in rotating fluids 
have focused on the inertial modes or wave attractors in 
closed containers ,'^^1'^ whereas less attention has been 
paid to localized inertial wave beams in effectively un- 
bounded systems. Inertial modes and attractors are gen- 
erated either from a disturbance of significant size com- 
pared to the container, or more classically from global 
forcing (precession or modulated angular velocity ).^^ ^ 
Localized inertial waves generated by a small distur- 
bance have been visualized from numerical simulations by 
Godeferd and Lollini,^^ and have been recently investi- 
gatedjusing Particle Image Velocimetry (PIV) by Messio 
et al}^ In this latter experiment, the geometrical prop- 
erties of the conical wavepacket emitted from a small 
oscillating disk was characterized, by means of velocity 
measurements restricted to a horizontal plane normal to 



the rotation axis, intersecting the wavepacket along an 
annulus. 

The weaker influence of rotation compared to stratifi- 
cation in most geophysical applications probably explains 
the limited number of references on inertial waves com- 
pared to the abundant literature on internal waves (see 
Voisin^^ and references therein). Another reason might 
be that quantitative laboratory experiments on rotating 
fluids are more delicate to perform than for stratified flu- 
ids. In particular, the velocity field of internal waves is 2 
components, whereas that of inertial waves is 3 compo- 
nents (although the geometry of the wave pattern may 
be two-dimensional). Moreover, only PIV is available 
for quantitative investigation of the wave structure for 
inertial waves, whereas other optical methods, such as 
shadowgraphy, or more recently synthetic Schlieren,!^ are 
also possible for internal waves. 

The purpose of this paper is to extend the results of 
Messio et al.^^^ using a newly designed rotating turntable, 
in which the velocity field can be measured over a large 
vertical field of view using a corotating PIV system. 
In the present experiment, the inertial wave is gener- 
ated by a thin cylindrical wavemaker, producing a two- 
dimensional cross-shaped wave beam, and special atten- 
tion is paid to the viscous spreading of the wave beam. 
The beam thickness and the vorticity decay are found 
to compare well with a similarity solution, analogous to 
the one derived by Thomas and Stevenson^ for internal 
waves. 



II. THEORETICAL BACKGROUND 
A. Geometry of the wave pattern 

A detailed description of the structure of a plane 
monochromatic inertial wave in an inviscid fluid can be 
found in Messio et al.^^^ and only the main properties 
are recalled here. We consider a fluid rotating at con- 
stant angular velocity = ^le^, where the direction ez 
of the reference frame (ex, ey, ez) is vertical (see Fig.jl]). 
Fluid particles forced to oscillate with a pulsation cj <2Vt 
describe anticyclonic circular trajectories in tilted planes. 
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FIG. 1. Geometry of an inertial wave beam emitted in an 
infinite medium from a localized oscillating cylindrical wave- 
maker invariant in the y-direction. 



A propagating wave defined by a wavevector k normal to 
these oscillating planes is solution of the linearized invis- 
cid equations, satisfying the dispersion relation: 



(J = 2Q.- k/k 21] cos ( 



In this relation, only the angle of k with respect to the 
rotation axis is prescribed, whereas its magnitude is set 
by the boundary conditions. For such anisotropic disper- 
sion relation, the phase velocity, c = ak/A;^, is normal to 
the group velocity,^ Cg = Vkcr (see Fig. IT]). 

If one now considers a wave forced by a thin horizon- 
tal velocity disturbance invariant in the Y direction, al- 
though the velocity field is still 3 components, the wave 
pattern is two-dimensional, varying only in the (X, Z) 
vertical plane. The wave pattern consists in four plane 
beams making angle ±6> with respect to the horizontal, 
drawing the famous St. Andrew cross familiar in the 
context of internal waves. ^ In the following, we consider 
only one of those four beams, with X > and Z > 0, 
and we define in Fig. [T] the associated local system of co- 
ordinates {ex^ey^Bz): The axis e^^ is in the direction of 
the group velocity, is directed along the wavevector k, 
and Gy = ey is along the wavemaker. 

Since the source is localized, a broad spectrum of 
wavevectors is excited, all aligned with e^. In an inviscid 
fluid, the interference of this infinite set of plane waves 
will cancel out everywhere except in the z = plan, re- 
sulting in a single thin oscillating sheet of fluid describing 
circular trajectory normal to e^. 



B. Viscous spreading 

In a viscous fluid, the energy of the wave beam is dissi- 
pated because of the shearing motion between oscillating 
planes. As the energy propagates away from the source, 
the larger wavenumbers will be damped first, so that the 
spectrum of the wave beam gradually concentrates to- 
wards lower wavenumbers, resulting in a spreading of the 
wave beam away from the source. 



Although the viscous attenuation of a single Fourier 
component yields a purely exponential decay, the atten- 
uation of a localized wave follows a power law with the 
distance from the source, which originates from the com- 
bined exponential attenuation of its Fourier components. 
A similarity solution for the viscous spreading of a wave 
beam has been derived by Thomas and Stevenson^ in the 
case of internal waves, and has been extended to the case 
of coupled internal-inertial waves by Peat.^ The deriva- 
tion in the case of a pure inertial wave is detailed in the 
Appendix, and we provide here only a qualitative argu- 
ment for the broadening of the wave beam. 

During a time t, the amplitude of a planar monochro- 
matic wave of wavevector k is damped by a factor 
e/e = exjp{—uk'^t) as it travels of a distance x = Cgt 
along the beam, where Cg is the group velocity. Using 
Cg = {2Q/k)sm6 = (cr/^) tan 6>, the attenuation factor 
writes 

Ck = exp(— -^^/c^x). 



(1) where we introduce the viscous lengthscale. 
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For a wave beam emitted from a thin linear source at x = 
0, an infinite set of plane waves are generated, and the 
energy of the largest wavenumbers will be preferentially 
attenuated as the wave propagates in the x direction. At 
a distance x from the source, the largest wavenumber, for 
which the energy has decayed by less than a given factor 
e*, is kmax = {i'^x)~^^^ Ine*. At this distance x, the wave 
beam thus results from the interference of the remaining 
plane waves of wavenumbers ranging from to kmax- Its 
thickness can be approximated by 5{x) ~ k^\^^ yielding 
d{x)/£ ~ {x/£)^^^. Mass conservation across a surface 
normal to the group velocity implies that the velocity 
amplitude of the wave must decrease as X 

More specifically, introducing the reduced transverse 
coordinate rj = zjx^l'^^^l'^ ^ a similarity solution exists 
for the velocity envelope. 
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where C/q is the velocity scale of the wave, and the analyt- 
ical expression of the non-dimensional envelope Eo{r]) is 
given in the Appendix. Similarly, the vorticity envelope 
can be written as 



uJo{x) W( 
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with Wq the vorticity scale. Although the normalized 
velocity envelope Eo{r])/EQ{0) has larger tails than the 
vorticity one Ei{r])/Ei{0), they turn out to be almost 
equal for r] < 4. The width at mid- height, defined such 
that Emir]ii2/2) £;^(0)/2, with m = 0, 1, is 7/1/2 - 



6.84 for both envelopes, so that the width of the beam in 
dimensional units is 
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Finite size effect of the source 



The similarity solution described here applies only in 
the case of a source of size much smaller than the viscous 
scale £. In the case of internal waves, Hurley and Kead}/^ 
(see also Flynn et al.^) have shown that for a source of 
large extent, vertically vibrated with a small amplitude, 
the wave could be approximately described as originating 
from two virtual sources, respectively located at the top 
and bottom of the disturbance. Following qualitatively 
this approach in the case of inertial waves forced by a 
horizontal cylinder of radius i?, the boundaries of the 



upper wave are given by 



R±5{x)/2, and that of 



the lower wave are given by z^^^^ = —R ± 5{x)/2. The 
lower boundary of the upper source intersects the upper 
boundary of the lower source at a distance Xi such that 



-up 



-down 



(xi), yielding 




Zo{t) = A cos{aot) 



Visualization 



FIG. 2. Schematic view of the experimental setup. The hori- 
zontal 8 mm diameter cylinder is oscillating vertically accord- 
ing to Zo{t) = A cos{ao t), with A — 2 mm and Go — 0.2 Hz. 
PIV measurements in a vertical plane (X^Z^ in the rotating 
frame are achieved by a vertical laser sheet and a camera at 
90°. 
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For large wavemakers {Rji ^ 0.025~^/^ 6.3), one 
has two distinct wave beams for x <C x^, and one sin- 
gle merged beam for x ^ xi. On the other hand, for 
smaller wavemakers, the merging of the two wave beams 
occur virtually inside the source, which can be effectively 
considered as a point source. In this case, the effective 
beam width far from the source may be simply written 
as 

b,^{x)-2R^b{x). (7) 



III. THE EXPERIMENT 
A. Experimental setup 

The experimental setup consists in a cubic glass tank, 
of sides 60 cm and filled with 54 cm of water (see Fig. [2]), 
mounted on the new precision rotating turntable "Gy- 
roflow", of 2 m in diameter. The angular velocity Vt of 
the turntable is set in the range 0.63 to 2.09 rad s~^, with 
relative fluctuations IS^tlVt less than 5 x 10~^. A cover is 
placed at the free surface, preventing from disturbances 
due to residual surface waves. The rotation of the fluid 
is set long before each experiment in order to avoid tran- 
sient spin-up recirculation flows and to achieve a clean 
solid body rotation. 

The wavemaker is a horizontal cylinder of radius R — 
4 mm and length L = 50 cm, hung at 33.5 cm below the 
cover by a thin vertical stem of 3 mm diameter. It is 



off-centered in order to increase the size of the studied 
wave beam in the quadrant X < and Z < 0. The ver- 
tical oscillation Zo(t) = A cos(<Jot), with A = 2 mm, is 
achieved by a step- motor, coupled to a circular camshaft 
which converts the rotation into a sinusoidal vertical os- 
cillation. In the present experiments, the wavemaker fre- 
quency is kept constant, equal to ao = 1.26 rad s~^, 
and the angular velocity of the turntable is used as the 
control parameter. This allows the velocity disturbance 
aoA = 2.5 mm s~^ to be fixed, whereas the angle of 
the inertial wave beam with respect to the horizontal, 
= cos~^(<Jo/2r^), is varied between and 72°. The 
velocity and vorticity profiles are examined at distances 
X between 30 and 300 mm from the wavemaker. The 
three-dimensional effects originating from the finite size 
L of the cylinder can be safely neglected since x < 0.6L. 
The Reynolds number based on the wavemaker velocity 
is Re = croA(2R) / V ~ 20, so that the flow in the vicin- 
ity of the wavemaker is essentially laminar. Except in 
Sec. IV B[ where the transient regime is described, mea- 
surements start after several wavemaker periods in order 
to achieve a steady state. 



For the forcing frequency do considered here, the 
boundary layer thickness is 8s = (y jcFo)^!^ 0.9 mm. 
This thickness also giv es the order of magnitude of the 
viscous length i = Ss/VtanO [see Eq. ([2])], for angles not 
too close to and 7r/2. The wavemaker radius being cho- 
sen such that R/i 4, the small source approximation is 



satisfied according to the criterion discussed in Sec. II C 
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B. PIV measurements 

Velocity fields in a vertical plane (X,Z) are measured 
using a 2D Particle Image Velocimetry system. The flow 
is seeded by 10 /im tracer particles, and illuminated by 
a vertical laser sheet, generated by a 140 mJ Nd:YAg 
pulsed laser. A vertical 43 x 43 cm^ field of view is ac- 
quired by a 2048 x 2048 pixels camera synchronized with 
the laser pulses. The field of view is set on the lower 
left wave beam. For each rotating rate, a set of 2000 im- 
ages is recorded, at a frequency of 2 Hz, representing 10 
images per wavemaker oscillation period. 

PIV computations are performed over successive im- 
ages, on 32 X 32 pixels interrogation windows with 50% 
overlap, leading to a spatial resolution of 3.4 mm.^^ In 
the following, the two quantities of interest are the ve- 
locity component obtained from the measured com- 
ponents ux and uz projected along the direction of the 
wave beam, and the vorticity component Uy^ normal to 
the measurement plane. 

The velocity along the wave beam typically decreases 
from 1 to 0.1 mm s~^, and is measured with a resolution 
of 0.02 mm s~^. Two sources of velocity noise are present, 
both of order of 0.2 mm s~^, originating from residual 
modulations of the angular velocity of the turntable, and 
from thermal convection effects due to a slight difference 
between the water and the room temperature. The resid- 
ual velocity modulations, of order of Lq AQ/2 (where Lq 
is the tank size and AQ 5 x 10"^^^) are readily re- 
moved by computing the phase- averaged velocity fields 
over the 200 periods of the wavemaker. Thermal con- 
vective motions, in the form of slowly drifting ascend- 
ing and descending columns, could be reduced but not 
completely suppressed by this phase-averaging, and rep- 
resent the main source of uncertainty in these experi- 
ments. However, the vorticity level associated to those 
convective motions appears to be negligible compared to 
the typical vorticity of the inert ial wave. Therefore, the 
vorticity profiles of the wave could be safely computed 
from the phase-averaged velocity fields. 
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FIG. 3. Close-up view of the phase averaged velocity (arrows) 
and vorticity Uy (color mapped) for an experiment performed 
at ao/2Q = 0.67. The black line shows the direction predicted 
by the dispersion relation cos^ = ao/2Q. 



(a) (P = 7t/5 (b) <?^=27r/5 
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IV. GENERAL PROPERTIES OF THE WAVE 
PATTERN 

A. Visualisation of the wave beams 

Figure [3] is a close-up view of the velocity and vortic- 
ity fields at cio/^Q = 0.67, showing the vorticity layers 
of alternating sign where the measured velocity is almost 
parallel to the ray direction e^^. The angle of the ray 
with respect to the horizontal accurately follows the pre- 
diction of the dispersion relation ([T]) , as illustrated by the 
black line. In Fig. |4ja) to (c), phase- averaged horizontal 
vorticity fields Uy are shown for three regularly sampled 
values of the phase. One can clearly see the location of 
the inertial wave inside a wavepacket that draws the clas- 
sical four-rays St. Andrew cross. The evolution of the 



FIG. 4. Phase averaged horizontal vorticity field Uy for 
ao/2Q = 0.67 at different phases: (a) = 7r/5, (b) = 27r/5 
and (c) (/) = 37r/5. The black line in (a) draws the direction 
predicted by the dispersion relation, (d) Vorticity envelope 
field (jjQ (see Sec.[v|. The dashed black and white lines shows 
the wave beam thickness predicted by the similarity solution 
[see Eq. 



vorticity field from Fig. [Ij^a) to (c) illustrates the prop- 
agation of the phase in directions normal to the rays. 
Some reflected wave beams of much smaller amplitude 
may also be distinguished on the background. 

Figure |5] compares the cosine of the measured angle 
with the theoretical value (Jo/2Q according to the pre- 
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(a) (7o/2n = 0.85 



(b) (Jo/2r2 = 0.80 



0.4 0.6 

FIG. 5. Cosine of the measured average beam angle, cos^, as 
a function of the frequency ratio cro/2Q. The hne shows the 
dispersion relation ([T]). Experimental uncertainties are of the 
order of the marker size. 



diction of the dispersion relation ([T]). The angle is deter- 
mined from the location of the maximum of the vorticity 
envelope (envelope computations described in Sec. V|), 
shown in Fig. [4]^d), averaged along the ray. An excellent 
agreement is obtained between the measurements and the 
theory, with a relative error less than 3%. 



B. Transient experiments 

In order to characterize the formation of the inertial 
wave pattern as the oscillation is started, a series of tran- 
sient experiments has been performed. In the case of a 
pure monochromatic plane wave, the front velocity of 
the wavepacket is simply given by the group velocity. 
However, in the case of a localized wave beam, since 
each Fourier component k travels with its own group 
velocity Cg = (<j//c)tan^, the shape of the wavepacket 
gradually evolves as the wave propagates. A rough esti- 
mate for the front velocity can be simply obtained from 
Vf ^ cr( A/27r) tan ^, where A is the apparent wavelength 
of the wave, simply estimated as twice the distance be- 
tween the locations of two successive vorticity extrema. 

Figure |6] shows spatiotemporal diagrams of the vortic- 
ity ujy{x, z = O^t) at the center of the beam as a function 
of the distance x from the wavemaker, for GqI^VL between 
0.85 and 0.50. Superimposed to these spatiotemporal 
images, we show the front velocity Vj o:^ cr (A/27r) tan ^, 
starting from x = at t = 0. Qualitative agreement 
with the spatiotemporal diagrams is obtained, indicat- 
ing that the propagation of the wave envelope is indeed 
compatible with this front velocity. 

Further quantitative estimate of the front velocity 
would require to extract the instantaneous wave enve- 
lope from those spatio-temporal diagrams, which is dif- 
ficult because the front velocity and the phase velocity 
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FIG. 6. Spatiotemporal representation of the vorticity ojy 
along the wave beam, where space is the distance x to the 
oscillating cylinder, for experiments performed at ao/2fl — 
0.85, 0.80, 0.75, 0.66, 0.60, 0.50. Black lines originating at 
(x = 0,t = 0) trace the front velocity Vf = a{\/ 2n) ta n^ es- 
timated from the apparent wavelength (see Sec. VA). t = 
corresponds to the start of the oscillation. 



are of the same order. This property actually prevents 
from a safe extraction of a longitudinal wavepacket enve- 
lope using standard temporal averaging over small time 
windows. 



C. Generation of harmonics 

Going back to steady waves, we now characterize the 
generation of higher order wave beams which take place 
at low forcing frequency. According to the dispersion re- 
lation, an harmonic wave of order n > 2 is allowed to 
develop whenever ncFo/2Vt < 1. Such harmonic waves 
of order n > 2 may originate either from a residual 
non-harmonic component of wavemaker oscillation pro- 
file Zo{t), or from inertial non-linear effects in the flow 
in the vicinity of the wavemaker, which can exist at the 
Reynolds number Re :^ 20 considered here. 

In Fig. [4j in which (Jo/2Q = 0.67, only the fundamen- 
tal wave (n = 1) can be seen. On the other hand, in 
Fig. [7|^a), in which ao/2Q = 0.43, a second harmonic 
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FIG. 7. (a) Phase averaged vorticity field ouy for an exper- 
iment performed at ao/2Q — 0.43, siiowing botii tiie fun- 
damental (n = 1) and the second harmonic (n = 2) wave 
beams. The corresponding frequency- filtered vorticity fields 
are extracted in (b) and (c). 



wave beam is clearly present, propagating at an angle 
closer to the horizontal, as expected from the dispersion 
relation. This is confirmed by Fig. [7|^b) and (c), showing 
the corresponding frequency-filtered phase averaged vor- 
ticity fields, in (b) for the fundamental n = 1 and in (c) 
for the second harmonics n = 2. 

In order to further characterize this generation of har- 
monics, we have performed a spectral analysis of the time 
series of the longitudinal velocity Ux(t)^ measured at a 
given distance xq — 100 mm from the source, at the 
center of each wave beam. The energy spectrum |i^crp, 
where Ua is the temporal Fourier transform of Ux(t\ is 
shown in Fig. [s] for the two cases (Jo/2Vt = 0.67 and 0.43. 
In both cases, the spectra are clearly dominated by the 
fundamental forcing frequency cTq. Two other peaks are 
also found, at a = Q and a = originating from the 
residual modulation of t he ang ular velocity of the plat- 
form, as discussed in Sec. Ill B| (the energy of those peaks 
is typically 3 to 10 times smaller than the fundamental 
one). As expected, no harmonic frequency nao {n > 2) 
is found for ao/2Q = 0.67 (see Fig. [S^i), but a second 
harmonic n = 2 is indeed present for cio/^Q = 0.43 (see 
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FIG. 8. Energy spectrum of the velocity time series measured 
at the center of the wave beam of interest, at a fixed distance 
xo — 100 mm from the wavemaker. (a) ao/2Q — 0.67, show- 
ing a single peak at the forcing frequency, (b) ao/2Q = 0.43, 
showing measurements performed in the fundamental beam 
n = 1 (in red/light gray) and in the second harmonic beam 
n = 2 (in blue/dark gray). In (a) and (b), insert shows the 
same spectrum in semi- logarithmic coordinates. Additional 
peaks are present at a/2Q = 0.5 and 1, originating from me- 
chanical noise of the rotating platform. 



Fig. [sJd). In this case, the energy ratio of the first to the 
second harmonics, each of them being measured at a dis- 
tance Xq = 100 mm from the source on the corresponding 
beam, is |'?i2crP/|'^crP — 0.036. As ao/2Q is further de- 
creased, the ratio |'?/2aP/|'^crP increases, reaching 0.05 for 
(Jo/^Q = 0.30, and even higher order harmonics emerge 
although with very weak amplitude. 



V. TEST OF THE SIMILARITY SOLUTION 

A. Velocity and vorticity envelopes 

We now focus on the dependence of the wavepacket 
shape and the viscous spreading of the wave beam with 
the distance x from the source. Figures [9ja) and (b) 
illustrate the shape of the phase- averaged velocity and 
vorticity profiles respectively, for two values of the phase 
(j)o and (/)o + 27r/5. The wavepacket envelopes are defined 
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FIG. 9. (a) Velocity envelope uo{xo, z) and two velocity pro- 
files Ux{xo, (j)) for two values of the phase 0, as a function of 
the transverse coordinate z at a fixed distance Xo — 100 mm 
from the wavemaker for cro/2Q = 0.67. (b) Corresponding 
vorticity envelope ujq{xq^ z) and vorticity profiles ujy{xQ^ z^cj)). 
5 is the envelope thickness at mid- height. In (a) and (b), the 
green/light gray curves correspond to the experimental en- 
velope, and the dotted curves correspond to the predictions 
of the similarity solutions normalized by the measured maxi- 
mum. 



Uq{x,z) = ^2{U:j,{x,Z,(j)Y)^ 

(and similarly for cjq), where is the average over all 
phases (j). Although the measured normalized envelopes 
compare well with the normalized envelopes predicted 
from the similarity solutions [£^^(7^)/£^^(0), with m = 
for the velocity and m = 1 for the vorticity], the agree- 
ment is actually better for the vorticity. This is prob- 
ably due to the velocity contamination originating from 
the residual angular velocity modulation of the platform 
and the slight thermal convection effects discussed in 
Sec. |IIIB[ The better defined vorticity envelopes actu- 
ally confirms that those velocity contaminations have a 
negligible vorticity contribution. For this reason, we will 
concentrate only on the vorticity field in the following. 
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FIG. 10. (a) Velocity envelope uq{xq^z) and three velocity 
profiles Ux{xq^ z^cj)) as a function of the transverse coordinate 
z at a fixed distance Xo = 70 mm from the wavemaker for 
ao/2Q = 1. (b) Corresponding vorticity envelope and vor- 
ticity profiles. The arrows indicate the time evolution of the 
profiles. The interference of the upward and downward wave 
beams produces a stationary wave pattern at z = with a 
velocity node and a vorticity maximum. 



It is worth to examine here the singular situation 
(Jo/^Q = 1, in which the similarity solution is no longer 
valid. In this situation, the phase velocity is strictly ver- 
tical and the group velocity vanishes. The upward and 
downward beams are expected to superimpose and gen- 
erate a stationary wave pattern in the horizontal plane 
Z = z 



0. Figure 10 



shows the velocity envelope 
uo{xo^ z) and three phase averaged profiles as a function 
of the transverse coordinate z. The observed wave is 
actually stationary at the center of the wavepacket (see 
the velocity node and vorticity maximum for z = 0), 
and shows outwards propagations on each side of the 
wavepacket. 

Returning to the standard situation cFo/2Vt < 1, the 
vorticity amplitude at a given location x is defined as 
the maximum of the vorticity envelope at the center of 
the beam, ujmax{x) = ujo{x,z = 0). The thickness of the 
wavepacket S{x) is defined from the width at mid- height 
of the envelope, such that 

{x)/2. 
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This beam thickness 5 depends both on the distance x 
from the source and on the viscous length £ [see Eqs. ([5| 
and In order to check those two dependencies, 5 is 
plotted in Fig. 
and in Fig. [IT 
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^a) as a function of x at fixed ao/2Q, 
Id) as a function of ao/2Q at fixed xq. 
The agreement with the effective wave beam thickness 
4fr = 2R^6Mi{x/iy^^ is correct, to within 10%, which 
justifies the simple analysis of merged beams originat- 
ing from the two virtual sources located at the top and 
bottom of the wavemaker. The oscillations of 6 proba- 
bly originate from the interaction of the principal wave 
beam with refiected ones. Figure [llja) also shows the 
apparent wavelength X{x) of the wave, simply defined as 
twice the distance between a maximum and a minimum 
of the phase- averaged vorticity profiles. This apparent 
wavelength turns out to be even closer to the expected 
lengthscale (5eff of Eq. ([7|), to within 4%, suggesting that 
A is less affected by the background noise than the beam 
thickness. A good agreement between both S and A and 
the prediction ^ is also obtained as ao/2Q (and hence £) 
is varied at fixed xq, as shown in Fig. [iTJb). Here again, 
the interaction with refiected wave beams is probably re- 
sponsible for the significant scatter in this figure. 



B. Decay of the vorticity envelope 

The decay of the vorticity amplitude ujmax{x) as a 
function of the distance x from the source is shown in 
Fig. 12 Taking the similarity solution Q at the center 
of the wave beam z = yields 
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FIG. 11. Wave beam thickness 6 (red o), and apparent wave- 
length A (blue □); (a), as a function of the distance x from the 
wavemaker for ao/2Q = 0.67; (b), as a function of ao/2Q at 
a distance Xo — 100 mm from the wavemaker. In both plots, 
the line shows the predicted effective wave beam thickness ^eff 

0. 



Letting the vorticity scale Wq as a free parameter, a 
power law is found to provide a good fit for the 

overall decay of ujmax{x). Some marked oscillations are 
however clearly visible, e.g. at x between 220 and 320 mm 
for ao/2n = 0.85 (blue □ in Fig.[l2|. Those oscillations 
appear at locations where reflected wave beams interact 
with the principal one, inducing modulations of the wave 
amplitude. This interpretation is confirmed by the fact 
that (i) the observed modulation has a wavelength of 
45 mm, which corresponds to the apparent wavelength 
of the wave, and that (ii) in Fig. [4j corresponding to 
Go 120. = 0.67, a modulation of the principal wave beam 
by a refiected can be clearly seen at a distance of about 
250 mm from the source. 

The vorticity scale Wq is theoretically related to 



rr* 
^0 



through the relation I^q 



the velocity scale 
{Ei{0)/Eo{0))U^/£ 0.506 /7o7^ (see the Appendix). 
Since the wavemaker velocity is ctqA, the velocity scale 
Uq is expected to write in the form aoAg{0), where the 
unknown function g{0) describes the forcing efficiency of 
the wavemaker. Accordingly, the forcing efficiency can 
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FIG. 12. Vorticity amplitude uJmax{x) as a function of the 
distance x from the wavemaker for ao/2Q = 0.67 (blue □) and 
ao/2Q = 0.85 (red o), and best fit with the law WS{x/£y 
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FIG. 13. Forcing efficiency g{0) defined from Eq. (|9| as a func- 
tion of ao/2Q. Squares and errorbars represents respectively 
the mean and the standard deviation for each ao/2Q ratio 
(reflecting the variability of uJmax along x). The black line 
shows a best fit according to Eq. (10), with go = 0.94 zb 0.10. 



be deduced from the vorticity data, by computing 



9(0) 



: (x/£)2/3 



0.506 aoA/£ 0.506 aoA/i 



(9) 



for each value of ao/2Q. Measurements of g{0) are plot- 
ted as a function of ao/2Q in Fig. [l3] As expected, this 
forcing efficiency decreases as GqI^VL is increased, i.e. as 
the wave beam becomes closer to the horizontal. In the 
limit (Jo 1 2^ 1, the vertically oscillating wavemaker be- 
comes indeed very inefficient to force the quasi-horizontal 
velocities of the wave. 

An analytical expression for the function g(d) would 
require to solve exactly the velocity field in the vicinity 
of the wavemaker, and in particular the coupling between 
the oscillating boundary layer and the wave far from the 
source, which is beyond the scope of this paper. In the 
case of a cylinder, a naive estimate of g(S) could however 
be obtained, assuming that the effective velocity forcing 
is simply given by the projection of the wavemaker ve- 
locity along the wave beam direction, yielding 



g{d) ^osin6> = go \ I 



\2QJ 



(10) 



with go a constant to be determined. A best fit of 
the experimental values of g{6) with this law leads to 
gQ :^ 0.94 ± 0.10 [see Fig. [l3], and reproduces well the 
decrease of g{0) as ao/2Q is increased. The fact that 
go is found close to 1 indicates that the inertial wave 
beam is essentially fed by the oscillating velocity field in 
the close vicinity of the wavemaker. The discrepancy at 
large forcing frequency may be due to the breakdown of 
the similarity solution as the angle approaches 0. 



VI. STOKES DRIFT 

We finally consider the possibility of a Stokes drift 
along the wavemaker which is expected because of the 
attenuation of the wave amplitude induced by the viscos- 
ity. A fluid particle in the inertial wave approximately 
describes a circular orbit. During this orbit, the parti- 
cle experiences a larger velocity along y when it is closer 
than when it is further from the wavemaker (see Fig. [T]), 
resulting in a net mass transport along y. This is simi- 
lar to the classical Stokes drift for water waves, which is 
horizontal because of the decay of the velocity magnitude 
with depth. Here the drift is expected in general in the 
direction given by x c^, since the viscous attenuation 
takes place along the group velocity Cg. 

Attempts to detect this effect have been carried out, 
from PIV measurements in vertical planes (Y^Z). Be- 
cause of the weakness of the considered drift, the mea- 
surements have been performed very close to the wave- 
maker, for X between 5 and 30 mm, where a stronger 
effect is expected. However, those attempts were not 
successful, probably because the drift, if present, is hid- 
den by the stronger fluid motions induced by the residual 
thermal convection columns, as discussed in Sec. IIIIB 



The magnitude of the expected Stokes drift cannot be 
easily inferred from the complex motion of the fluid par- 
ticles close to the wavemaker. An estimate could however 
be obtained in the far field, from the similarity solution 
of the wave beam. We consider for simplicity a particle 
lying at the center of the wave beam (z = 0), at a mean 
distance xo from the source, describing approximate cir- 
cles of gyration radius a |u(xo)|/cr in the tilted plan 
(x^y). The expected drift velocity can be approxi- 
mated by computing the velocity difference between the 
two extreme points xo — a and xo+a of the orbit, yielding, 
to first order in a/xo to^ 



VSy{xo) - Y 



ax, 



5/3 



(11) 



The steep decrease as x^^^^ confirms that the drift should 
be essentially present close to the wavemaker. Although 
this formula is expected to apply only in the far-field 
wave (typically for xo > 40^, see the Appendix), its ex- 
trapolation close to the wavemaker, for xo — 10£ o:^ 2R, 
gives vsy ^0.1 mm s~^. This expected drift velocity 
is about 10% of the wave velocity at the same location, 
but it turns out to remain smaller than the velocity con- 
tamination due to the thermal convection columns. Al- 
though the phase averaging proved to be efficient to ex- 
tract the inertial wave field from the measured velocity 
field - thanks to a frequency separation between convec- 
tion effects and inertial wave -, it fails here to extract 
the much weaker velocity signal expected from this drift, 
since it is zero frequency and hence mixed with the very 
low frequency convective noise. 
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VII. CONCLUSION 
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In this paper, Particle Image Velocimetry measure- 
ments have been used to provide quantitative insight into 
the structure of the inertial wave emitted by a verticahy 
oscihating horizontal cylinder in a rotating fluid. Large 
vertical fields of view could be achieved thanks to a new 
rotating platform, allowing for direct visualization of the 
cross-shaped St. Andrew wave pattern. 

It must be noted that performing accurate PIV mea- 
surements of the very weak signal of an inertial wave is a 
challenging task. In spite of the high stability of the an- 
gular velocity of the platform {Aft /ft < 5 10~^), the ve- 
locity signal-to- noise remains moderate here. Addition- 
ally, slowly drifting vertical columns are present because 
of residual thermal convection effects, and are found to 
account for most of velocity noise in these experiments. 
Those thermal convection effects are very difficult to 
avoid in large containers, even in an approximately ther- 
malized room. However, this noise can be significantly 
reduced by a phase-averaging over a large number of os- 
cillation periods. This concern is not present for inter- 
nal waves in stratified fluids, because residual thermal 
motions are inhibited by the stable stratification. This 
emphasizes the intrinsic difficulty of experimental inves- 
tigation of inertial waves, in contrast to internal waves 
which have been the subject of a number of studies. 

In this article, emphasis has been given on the spread- 
ing of the inertial wave beam induced by viscous dis- 
sipation. The attenuation of a two-dimensional wave 
beam emitted from a linear source is purely viscous, 
whereas it combines viscous and geometrical effects in 
the case of a conical wave emitted from a point source. 
The linear theory presented in this paper is derived un- 
der the classical boundary-layer assumptions first intro- 
duced by Thomas and Stevenson^ for two-dimensional 
internal waves in stratified fluids. The measured thick- 
ening of the wave beam and the decay of the vorticity 
envelope are quantitatively fitted by the scaling laws of 
the similarity solutions of this linear theory, S{x) ~ x^^^ 
and 0Jmax{x) ~ where x is the distance from the 

source. More precisely, we have shown that the ampli- 
tude of the vorticity envelope could be correctly predicted 
from the velocity disturbance induced by the wavemaker, 
by introducing a simple forcing efficiency function g{6)^ 
where is the angle of the wave beam. 

Finally, it is shown that an attenuated inertial wave 
beam should in principle generate a Stokes drift along 
the wavemaker, in the direction given by 11 x c^, where 
Cg is the group velocity. However, in spite of the high 
precision of the rotating platform and the PIV measure- 
ments, attempts to detect this drift were not successful 
in the present configuration. Velocity fluctuations in- 
duced by thermal convection effects probably hide this 
slight mean drift velocity, suggesting that an improved 
experiment with a very carefully controlled temperature 
stability would be necessary to detect this very weak ef- 
fect. 



We acknowledge A. Aubertin, L. Auffray, C. Borget, 
G.-J. Michon and R. Pidoux for experimental help, and 
T. Dauxois, L. Gostiaux, M. Mercier, C. Morize, M. 
Rabaud and B. Voisin for fruitful discussions. The new 
rotating platform "Gyroflow" was funded by the ANR 
(grant no. 06-BLAN-0363-01 "HiSpeedPIV"), and the 
"Triangle de la Physique" . 



Appendix A: Similarity solution for a viscous planar 
inertial wave 



In this appendix, we derive the similarity solution for a 
viscous planar inertial wave, following the procedure first 
described by Thomas and StevensorP for internal waves. 

We consider the inertial wave emitted from a thin lin- 
ear disturbance invariant along the Y axis and oscillating 
along Z with a pulsation <j in a viscous fluid rotating at 
angular velocity ft = Qez. Since the linear source is in- 
variant along y, so will the wave beams, and the energy 
propagates in the (X, Z) plan. In the following, we con- 
sider only the wave beam propagating along X > and 
Z >0. 

The linearized vorticity equation is 

dtuj = (217 • V)u + i/V^cj. 

Recasting the problem in the tilted frame of the wave, 
(ex^ey^ez)^ with = ey, and e^^ tilted of an angle 
= cos~^{a/2Q) with the horizontal, one has ft = 
Q (sin Gx -\- cos e^), so that {2ft • V) = 2(7 (sinO dx + 
cosOdz) = a{tdin6dx + dz). Assuming that the 
flow inside the wave beam is quasi-parallel (boundary 



layer approximation), i.e., such that 



> Uz 



\^x\A^y\ ^ \^z\^ and \/'^ c:^ d'^^ the linearized vorticity 
equation reduces to 



dtUJx = cr{tdinOdx + dz)ux + i^d^uJx, 
dtOJy = cr{tdinOdx + dz)uy + ud'^ujy. 



(Al) 
(A2) 



We introduce the complex velocity and vorticity fields in 
the (x, y) plan as 



U = Ux -\- iu. 



w 



■ lUJ, 



Since, within the quasi-parallel approximation, one has 
W = idzU^ the combination (Al)+z(A2) yields 



idtdzU = a{t8inOdx + dz)U + iud^U. (A3) 



Searching solutions in the form U = U^e Eq. (A3) 
becomes 



(A4) 



whe re w e have introduced the viscous scale I ([2|. 



variable: 



Eq. (|A4|) admits similarity solutions as a function of the 
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which are of the form: 



Uo{x,z) = Uo[-] fir,), 



1/3 



(A6) 



and where Uq is a velocity scale and f{r]) a non- 
dimensional complex function of the reduced transverse 
coordinate rj. Plugging such similarity solution (A6) into 



Eq. (A4) shows that /(t^) is solution of the ordinary dif- 



ferential equation 



3/"' + i(/ + r//')=0, 



(A7) 



which is identical to the equation (16) derived by Thomas 
and Stevenson'' for the pressure field of internal waves. 
Following their development, we introduce the family of 
functions defined through 



fmiji) — H~ ^<^r) 



Jo 



(A8) 



where Cm a nd S m are real, and such that fo{v) is a solu- 
tion of Eq. ( [A7| ). 

The velocity in the plan of the wave beam is therefore 
given hy Ux = ^{U} and Uy = ^{U}, leading to 



Eo{0) \X 



1/3 



[co{r]) cos(crt) + so{r]) sin(crt)] 



^o(O) 



1/3 



[soiv) cos(crt) - co{r]) sin(crt)] 



with Uq = Eo{0)Uo 0.893^/o where we introduce the 
family of envelopes Emiv) = \fmiv)\ = (c^ + ^m)^^^ fo^" 
m = 0, 1. 

Similarly, the vorticity in the plan of the wave beam is 



^{W} and 



^{VI^}, so that 



Ei{0) \x 



Ei{0) \x 



2/3 



2/3 



-ci{r]) cos(crt) — si{r]) sin(crt)] , 



-si{r]) cos((jt) + Ci{r]) sm{(jt)] . 



with Wo* = [Ei{0)/Eo{0)] US/i - 0.506/7o7^. 

The velocity and vorticity envelopes, defined as uq = 
{«) + {uDf^ and CO = {{i^D + K))V2, where (•) is 
the time-average over one wave period, are given by 



Uq = Uq 



1/3 



Eoiv) 



XI Eq{0)' 



The two normalized envelopes Emiv) / ^m{^) are com- 
pared in Fig. 14 Interestingly, they closely coincide up 



to 77 c± 4, but the vorticity envelope decreases much more 
rapidly than the velocity envelope as 7^ ^ 00 (one has 
Em oc 1/7/^+^ for 7/ > 1). 
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FIG. 14. Normalized velocity ( — , m = 0) and vorticity (- -, 
m = 1) envelopes of the similarity solutions. 



It is interesting to note that velocity and vorticity in 
the present analysis are analogous to the pressure and 
velocity in the analysis of Thomas and Stevenson.^ One 
consequence is that the lateral decay of the velocity en- 
velope is sharper for an internal wave (as l/r/^) than for 
an inertial wave (as l/rj). 

Finally, the z component of the velocity is obtained 
using incompressibility {dxUx + dzUz = 0), yielding 



1 



3 \x 



2/3 



1 z 
3 X 



(A9) 



This result shows that the streamlines projected in the 
(x,z) plan are along lines of constant phase (i.e., con- 
stant 7^). As a consequence, a particle trajectory is an 
approximate circle, projected on curved surfaces, invari- 
ant along 7/, and such that z = rf P'l'^x^l'^ ^ with rf given 
by the initial location. 

The thickness 7^1/2 of the velocity and vorticity en- 
velopes are defined such that £^^(7/1/2/2) = £^rn(0)/2. 
Those thicknesses turn out to be almost equal: 7/1/2 — 
6.841 for 7n = and 771/2 — 6.834 for 7n = 1. In dimen- 
sional units, the wave thickness is thus given by Eq. ([s]). 
Interestingly, the envelope of is given by r]Eo{r])^ which 
tends towards 1 as 77 ^ 00, so that no thickness could be 
defined for 77^. 

Finally, we note that the quasi-parallel approximation 
used in the present analysis is satisfied for Itx^I/I'^^cI ^ 
1. Using Itx^I/I'^xI = V^^^ {i/x)^^^ /3, and evaluating the 
envelope ratio at the boundary of the wave, i.e., for r] = 
7/1/2/2 3.42, this criterion is satisfied within 10% for 
X > 38£. 
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